Measurement of scaling laws for shock waves in thermal nonlocal media 
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We are able to detect the details of spatial optical collisionless wave-breaking through the high aperture 
imaging of a beam suffering shock in a fluorescent nonlinear nonlocal thermal medium. This allows us to 
directly measure how nonlocality and nonlinearity affect the point of shock formation and compare results with 
numerical simulations. © 2013 Optical Society of America 
OCIS codes: 000.0000, 999.9999. 



Shock waves (SW) are characterized by an abrupt wave- 
profile, corresponding to discontinuous solutions of hy- 
perbolic partial differential equations [1]. Dissipative SW 
display a smooth wave front, due to viscous damp- 
ing, whereas dispersive, or collisionless, SW exhibit a 
characteristic oscillatory front. The latter are expected 
in Hamiltonian universal models for nonlinear media, 
such as the nonlinear Schrodinger (NLS) equation, in 
the weakly dispersive regime where hydrodynamical ap- 
proximations hold true [2-4]. Dispersion regularizes the 
shock, determining the onset of oscillations that appear 
near the wavebreaking point, and expand from there [2], 
as investigated in ion- acoustic waves [5], in fibers [6,7], 
in diffracting optical beams [8] or in Bose-Einstein con- 
densates [9]. Nonlocal phenomena [10-14] substantially 
alter dispersive SW [15-17]. For propagation in highly 
nonlocal, strongly defocusing, thermal media, previous 
work has been based on approximating the 3D dimen- 
sional Fourier equation for the temperature profile by 
a nonlocal nonlinearity with a finite degree of nonlocal- 
ity [16-18]. 

In this Letter we report on the quantitative investiga- 
tion of the formation of dispersive shocks from standard 
laser beams (Gaussian TEMqo) in thermal media. The 
key to our experiments is the direct observation of the 
intensity profile of the shock wave in its formation and 
evolution through the use of a fluorescent medium and 
a high aperture microscope. We quantify the shock po- 
sition along propagation and show how it depends on 
nonlinearity and nonlocality. Our results not only con- 
firm the finite degree of nonlocality of thermal media but 
actually provide a strategy for its direct measurement. 

In experiments we use aqueous solutions of Rho- 
damineB (RhB) displaying a thermal defocusing effect 
due to light absorption. A continuous-wave solid state 
laser at wavelength A=532 nm is focused on the in- 
put facet of the sample at the beam radius waist wq. 
We work at two different wq (18/im and and 25/im) to 
change both nonlinearity and nonlocality and two differ- 
ent RhB concentrations (0.1 mM, sample A; 0.067 mM, 



sample B) that vary the nonlinearity. For the visual- 
ization of light along its propagation, samples are put 
in 1cm X 1cm X 3cm glass cells (propagation along 1cm), 
and top images are collected through a microscope and 
recorded by a 1024 x 1392 pixel CCD camera. Conversely, 
to detect light transmitted at the exit of the samples, we 
use 1mm X 1cm X 3cm glass cells with propagation along 
the 1mm vertical direction (parallel to gravity) so as to 
moderate the effect of heat convection. Transverse im- 
ages of the beam intensity distributions at the exit facet 
of the cell are collected through a lens and recorded by 
the CCD camera. Images from the two experimental con- 
figurations are reproduced with comparable length-scale. 

Numerical simulations are also performed to verify 
the experimental measurements and to clarify the role 
of nonlocality. By using a split-step beam propagation 
method, we numerically integrate the following nonlocal 
NLS equation (details in [16,17]): 
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where ip = ^/x^d^ is the normalized field envelope 
(/ = \E\'^ is the intensity), 6 = koAnLd^ with An the 
nonlinear refractive index, Ld = koUoWQ the diffraction 
length, X = ^0^2, ko = 27r/A, and no the refractive in- 
dex, a = aoLd is the normalized loss coefficient defined 
in terms of the power absorption coefficient ao. Trans- 
verse spatial coordinates are scaled diS (,,r] = x/wo,y /wq, 
while C = z/Ld. Here a is the degree of nonlocality, used 
as a fitting parameter below. Figures la-c show three col- 
lected images of the transverse distribution of the beam 
intensity along z, for wq = 25/im and three beam powers 
P in sample A. In Fig. la, corresponding to the lowest 
P, the beam shape does not change (diffraction is weak) 
along propagation and the shock does not occur. The 
presence of oscillatory SW can be observed in the other 
two images reported in Fig. lb and Fig. Ic at higher 
powers, and in Figs, le-g where the 2-D profiles of the 
transmitted beam (on x-y plane) at the exit facet are 
reported. The x-profiles display post-shock rings that 
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increase in number and visibility with P, with outer 
rings being more intense than the inner ones, as typ- 
ical for dispersive SW from Gaussian beams in nonlo- 
cal media. Figures Id and Ih present beam propagation 
and transmission results obtained from simulations, per- 
formed by using input parameters identical to those in 
experiments and a = 0.2. Shock occurs in the regime 
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Fig. 1. (Color online) (a-d) Beam propagation as ob- 
served from top fluorescence emission for three different 
input powers (a-c) and calculated by numerical simula- 
tion (d); (e-g) corresponding transmitted intensity in the 
experiments and (h) in the simulations. Superimposed 
curves give the intensity x-profile. 

Lni <C La^Ld [2-4], with = l/o^o the absorption 
length, and Lni = l/xH^^^)peak is the nonlinear length 
(at the highest powers in our experiment, Lni/La ~ 
0.002 and, in terms of the smallness parameter defined 
in Refs. [16,17], e = ^jLni/Ld - 0.015). This means 
that the involved length scales are such that the beam is 
mainly affected by the nonlinear defocusing effect rather 
than by absorption and diffraction. In this regime, as 
a first approximation, the beam phase depends on the 
intensity profile as 0(^, C) = ^^'^{H^^v))- dis- 
cussed in Ref. [16], SW originates in correspondence with 
the chirp singularity \d(f)/d(,\ oo from which the shock 
distance can be analytically determined in the hydrody- 
namical approximation. For a generic degree of nonlo- 
cality and by assuming, for the sake of simplicity, a one- 
dimensional generic response function i^(^ — ^') yielding 
the index change as An = / d^' K{^-^')I{^' , r]) [Eqs. (1) 
imply a 2-D Lorentzian kernel K as one can easily verify 
by taking the Fourier transform of the second equation] , 
we have: 
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and the maximum of grad(0) along z (or () gives the 
shock position. Since K is ^^-independent the shock- 
point Zs = C^d can be directly estimated both in ex- 
periments and simulations as the point of largest steep- 
ness of beam intensity, as detailed below. We note that 
in the highly nonlocal regime, K{£^) = Kq^ and corre- 
spondingly d^(j) = {ko(Ko/no) L d^I = 0, i.e., as the de- 



gree of nonlocality increases, the response function filters 
out all the edge frequencies in the broadened intensity 
spectrum in proximity of the shock-point, thus delaying 
it. In the highly nonlocal limit, the mean- field approx- 
imation fully smooths out all the singularities and the 
shock disappears. To estimate the point of shock for- 
mation both in experiments and simulations, we calcu- 
late the X— derivative dxln of the intensity normalized 
to its peak value, say In{x^z). We define the steepness 
S{z) as the maximum over x (at any z) of the deriva- 
tive S{z) = mdiXx[dxlN{x^ z)]. Then the shock distance 
Zs is defined as the distance z corresponding to the max- 
imum steepness, [mdi-Xz{S{z)) z = Zs]. We follow this 
procedure for the whole set of images. In Fig. 2 we show 
representative profiles along with their derivatives at the 
input (Fig. 2a), right at the shock position (Fig. 2b), and 
beyond Zs (Fig. 2c), the steepness resulting the highest 
in Fig. 2b where the shock occurs. The behavior of the 
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Fig. 2. (Color online) Normalized intensity profiles In 
(black line) along with derivative {OxIn^ red line) for 
three different z at P=600mW: (a) input, z = 0; (b) 
z = Zs = 273/im, shock point; (c) z = 850/im, post- 
shock. 

steepness S along z as obtained from experimental data 
is illustrated in Fig. 3 for three power values. The shock 
distance Zg is given by the position of the peak of the 
ten-degree polynomial fitting curve (solid curves through 
data). In Fig. 3 it is clear that this distance is shorter 
for increasing power, meaning that larger laser intensi- 
ties speed up the shock formation process. We report in 
Fig. 4 the measured shock distance versus laser power, 
comparing with simulations. In particular Fig. 4a com- 
pares the shock point behavior for two beam waists in the 
same sample (B), whereas Fig. 4b compares the behavior 
in the two samples at a given beam waist (25/im). These 
data show quantitatively how both increasing input in- 
tensity or dye concentration speeds up shock formation 
because both result in a larger nonlinear response. In 
the simulations, we make use of the measured values of 
the physical parameters: no = 1.3, ao= 900 (A) and 
600 (B) m-^ a^L= 1.3x10-^ (A) and 1.2x10"^ (B) 
m W-^ ; n2 = 4 X 10"^^ (A) and 2.7 x 10"^^ (B) 
m^ W~^, while a is used as a free fitting parameter. 
The best-fit of numerical simulations with the experi- 



2 



(A 



J2 



250mW ■ 



460mW 




400 800 

z(nm) 



1200 



Fig. 3. (Color online) Steepness S{z) of transverse pro- 
files as function of z for three different powers. The point 
of maximum corresponds to shock distance Zg. 



mental data in Fig. 4a, yields a degree of nonlocality 
a = 0.63 {wq = 18/im) and 0.3 {wq = 25/im). Following 
Ref. [16], the nonlocality degree should scale with the in- 
verse of beam waist. The observed slight discrepancy is 
attributed to the theoretical two-dimensional heat equa- 
tion used to model the actual temperature profile in our 
sample. It is known from the hydrodynamic limit that 
the shock distance scales with power P according to the 
law Zs oc P~^, with 7 = 0.5 [16, 17]. This power law is 
confirmed by our experiments and simulations of Eqs. 
(1), though with deviations in the exponent. The exper- 
imental and numerical data given in Fig. 4, once fitted 
with 7 as a free parameter, yield the following values: on 
sample A, 7 = 0.88 {wq = 25/im), 7 = 1.0 {wq = 18/im) 
and on sample B, 7 = 0.70 {wq = 25/im), 7 = 0.72 
{wq = 18/im). The values of 7 from the simulations are: 
on sample A, 7 = 0.43 {wq = 25/im), on sample B: 
7 = 0.57 {wq = 25/im) and 7 = 0.43 {wq = 18/im). 
The discrepancy between simulations and experiments 
is, once again, ascribed to the limits of validity of the 2D 
reduction of the heat equation adopted in Eqs. (1). 

In conclusion, we use fluorescence imaging to directly 
monitor spatial wave-breaking. The full intensity dis- 
tribution map allows us to experimentally determine the 
shock position for different nonlinearity and nonlocality. 
The comparison with numerical simulations of the non- 
local NLS model enables to quantify the limits of validity 
of the 2D reduction of the heat equation. 
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Fig. 4. (Color online) (a) Shock position Zs versus in- 
put power for sample B at two different beam waists: 

■ 25/im , • 18/im; (b) Shock position vs. input power 
for beam waist=25/im and different dye concentrations: 

■ 0.067 mM, #0.1 mM. Experiments are denoted by 
black (dark) and simulations by red (clear) symbols. 
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